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1 Introduction 



The production of quarkonium states in high energy processes has recently attracted a lot 
of theoretical and experimental interest. The detection of J/0's plays a fundamental role 
in the study of B physics, because some of the most interesting final states of B decays do 
contain a J/ip. It is therefore important to have a good understanding of all other possible 
sources of J/0's, in particular the direct production. Furthermore, the large production 
cross sections and the relative ease with which J/0's can be triggered on even at small 
values of transverse momentum, make their observation a powerful tool to study hard 
phenomena in regions of small x, which are otherwise inaccessible with the standard hard 
probes (jets and vector bosons) used in high energy hadronic collisions. Since we expect 
J/0's to be mostly produced via gluon-gluon fusion, if a solid theory existed it could be 
used to get the best and most direct measurements of the gluon density of the proton at 
small x. The data are there, plentiful! 

Production models have existed for several years (see ref. U for a comprehensive review 
and references). However, it is only with the advent of the wealth of data from the high 
energy hadronic colliders [2-6] that significant tests of the theory have become possible, 
thanks to the big lever arm in CM energy relative to the fixed-target experiments, and 
thanks to the wide range in transverse momenta that can now be probed. The comparison 
of these data with the models available up to a couple of years ago has shown dramatic 
discrepancies, the most striking one (theory predicting a factor of 50 fewer prompt ip' 
than measured by CDF j|, [7|) having become known as the "CDF anomaly". 

Attempts to explain the features of these data have led in the past couple of years to 
a much deeper theoretical understanding of the mechanisms of quarkonium production. 
None of these developments would have been possible without the significant achievements 
of the experiments, to which most of the credit should go. In this review I summarize the 
evolution of the theoretical models towards what we can consider today as the seed of a 
theory of quarkonium production based on QCD. I will not have time to cover, however, 
the series of papers [EJ that proposed the existence of new exotic charmonium states to 
resolve the conflict between data and theory. 

The language I use is inspired by a space-time picture of the production process, and in 
my view should appeal to the generic reader for its simplicity. It should be kept in mind, 
however, that most of the qualitative statements that will be made can be rephrased in 
more rigorous terms, as discussed in the references that will be quoted. 

2 Quarkonium Production for Pedestrians 

Production of quarkonium represents a challenging theoretical problem, as it requires 
some understanding of the non-perturbative dynamics responsible for the formation of 
the bound state. The problem could be made easier by assuming the existence of some 
factorization theorem that allows the separation of the dynamics of the production of the 
heavy quark pair from its evolution into a bound state. The reason why this assumption 
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would make things simpler is that we could, in the presence of factorization, parametrize 
the non-perturbative part in a universal fashion in terms of a limited set of parameters: 
having been determined once (for example by fitting some set of data), these can then 
be used to perform predictions. The assumption of factorization is reasonable, since the 
time scales associated to the two phenomena are significantly different: in the production 
of the quark pair, the relevant time scale is the inverse of the mass of the heavy quark, 
or of its transverse momentum in the case of high-pr production. In the formation of the 
bound state, the important time scale is of the order of the inverse of the quarkonium 
binding energy, i.e. something of the order of 1/Aqcd- Therefore, by the time the bound 
state starts forming, the memory of what the source of the heavy quark pair was has been 
lost. However, the quarkonium state has well defined quantum numbers, and one might 
suspect that only heavy quark pairs prepared by the hard process in specific states have 
a chance to eventually evolve into a given bound state. Selection rules could therefore 
prevent the loss of memory, and spoil factorization. 

Determining to which extent, and in which precise form, factorization holds, is there- 
fore the primary challenge that we are faced with when formulating a production model. 
I will now shortly describe the two most popular models which have been proposed in 
the past, the so-called "colour evaporation" (CEM, ||) and "colour singlet" (CSM, fllPp 
models. They are based on orthogonal assumptions about the validity of factorization 
and, needless to say, lead to significantly different predictions. In a later section, I will 
show how these models evolved in the recent past into a more sophisticated approach, de- 
veloped within QCD, which can apparently explain the main features of the data currently 
available. 

2.1 Colour Evaporation 

In the colour evaporation (also known as local-duality) approach, factorization is assumed 
to hold strictly. Differential distributions for the production of a given quarkonium state 
H are assumed to be proportional to the production rate of a pair of heavy quarks with 
invariant mass in the range Ivhq < vhqq < 2m d'- 



Here Q is the heavy quark that forms the bound state, uiq is its mass and mu is the mass 
of the lightest meson carrying open flavour Q, i.e. the D meson in the case of charm, or 
the B meson in the case of bottom; rriQQ is the invariant mass of the produced heavy quark 
pair. The justification of this model stems from the assumption that only quark pairs 
below the threshold for production of open flavour can possibly bind into a quarkonium 
state, and that provided the quark pair has mass below this threshold, the correct quantum 
numbers for H will be recovered via non-perturbative emission {evaporation) of very soft 
gluons. The constant Ah, with Ah < 1, depends on the state H we are interested in, but 
is otherwise independent of the pt of H, and of the nature of beam and target. Therefore, 
while the model cannot estimate absolute cross sections, it however predicts their pt and 
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y/s dependence and it predicts ratios of production rates and distributions of different 
states to be a constant. 

As an example, consider production at large px' in this case the dominant mechanism 
for the production of heavy quark pairs with invariant mass close to threshold is the 
splitting of a high-p^ gluon: g — > QQ. The probability for this splitting to take place can 
be calculated in leading order (LO) QCD as: 



dProb a s 



dm 2 QQ 6tt m 2 QQ ' 



(2) 



and therefore: 



a s i 

dp? J * H hniQ Q ® \ dpff J Qir m 2 



QQ 



A H — ( da ( 9 A i g I mD \ = A H — ( da ( 9 A — (3) 
H 3tt \ dp% J \m Q J H 3n \ dp% J m H ' 

where we defined e = rri£, — raq <C ran- 

Attempts can be made to estimate what the relative values of Ah for different states 
H should be. For example, one could naively assume that a quark pair will evolve with a 
fixed probability into the closest state H with mass tuh < ^qq- In this case, 

da(H) oc (2J+1)— , (4) 

where e# is the mass splitting between adjacent states. In the case of charmonium, we 
would have for example m x — ~ 300 MeV, m^i — m x ~ 200 MeV, 2m D — tum ~ 
100 MeV, so that: 

*W) 



<t(Xj) 2J+1 



1/3 (5) 

(6) 



a{ X J>) 2J'+1 

J2 °{Xj) x BR( X j - VOM^O ~ 0.3. (7) 
j 

It is interesting to compare these naive predictions with data. From the publications of 
E705 ||rTJ (300 GeV beams of pions or protons on nuclei), the following ratios can be 
extracted: 

~~ 77Y ~ 0.25 8 



£ a( X j) x BR( X j - ^)/(r(ip) ~ 0.6 (9) 
j 
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The first result is also consistent with data from E789 [|12|], if the fraction of J/0's from \ 
decays is assumed to be of the order of 0.3. 

Results presented by DO at this Conference indicate a fraction of J/ip's from x's of 
approximately 0.4, however with a strong p? dependence. Likewise, CDF reported a 
value for this fraction of approximately 0.3, again with some pt dependence. CDF also 
measured the fraction of xi production relative to the total of \i + X2, of the order of 0.5. 
All of these numbers come quite close to the naive expectations given in eqs. (H)-©. The 
measurement of the total cross sections for T production by CDF |1 has also been shown 
to be consistent with the a/s dependence predicted by the colour evaporation model |13|] . 

Therefore, while not completely satisfactory and not a real theory, it is clear that the 
colour evaporation model presents some features of universality that are consistent with 
the observed data, and that should therefore be maintained by the final theory, however 
complicated this might be. 



2.2 The Colour Singlet Model 

The colour singlet model (TIJ, at least in its early formulation, emphasizes more the 
constraints imposed by the colour and spin selection rules. In the CSM, one projects the 
amplitude for the production of a heavy quark pair directly onto a state which has the 
right quantum numbers to form a given quarkonium state. This projection singles out 
only the right combinations of colour and spin required, and allows the absorption all the 
non-perturbative physics of the confinement into a single parameter, namely the value of 
the wave function (or derivatives thereof) of the quarkonium state at the origin. 

For example, in the case of 3 Si production this can be achieved by evaluating the 



following expression ]10|: 



Mty(P)) = ^£L^6,(P)Tr[^V(P + m)] (10) 

where M(ip) is the matrix element for the production of a 3 Si state of momentum P and 
mass m, is its polarization vector, R(0) is the wave function at the origin, and O l i is 
the matrix element for the production of the heavy flavour pair (i and j being the colour 
indices of quark and antiquark), with the constraint that the relative momentum between 
the two quarks be 0. Similar projection operators can be evaluated for any 2S+1 Lj state 
I(J, and convoluted with heavy quark production matrix elements calculable in QCD. 



Absolute predictions can therefore be made for the production rates, once we introduce 
the values of the wave functions that can be extracted from potential models or directly 
from the data on quarkonium decay widths. 

In the case of hadronic collisions, and at the leading order in a s , namely a 3 , it is easy 
to show that the only diagrams that are relevant for the production of a 3 Si state are 
those shown in Fig |T]a. At large transverse momentum, a simple kinematical analysis of 
the momentum flow in the two internal quark propagators shows that they are off-shell 
by approximately q 2 = — (4rriQ 2 +p^). As a result, the square of the matrix element 
will behave at large pt like l/pj<- This is a much steeper fall than that of standard hard 
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Figure 1: Sample graphs contributing to quarkonium production in the colour 
singlet model, (a): Leading order 3 Si. (b): Leading order 3 Pj. (c): Fragmenta- 
tion graph for 3 Si. The propagators crossed with a star are those with q 2 ~ — pjr 



processes, such as jet production, where the behaviour is that typical of a gluon or quark 
exchange, i.e. 1/pt- A similar analysis can be done in the case of large-pr production 
of x states, which is dominated by diagrams like the one shown in Fig |l]b. Here helicity 
conservation at the triple gluon vertex causes the internal gluon propagator to behave like 
1/Pt, and the total amplitude squared is therefore only suppressed by a factor of l/p?- 
Therefore the pr distributions of J/ip and \ as predicted at the LO by the CSM are totally 
different, contrary to what is assumed in the CEM. 

Such a steep px dependence has been proved to be inconsistent with the Tevatron 
data, where the accessible range of px is very large Aside from the technical details 
of how this behaviour arises from the diagrams, there is a simple reason why the CSM 
predicts such a strong suppression of high-pr quarkonium production at LO. In the CSM, 
one forces the heavy quark pair to be in the right state already at time scales significantly 
earlier than the time at which the formation of the bound state starts. This requirement 
produces a strong penalization in rate, which becomes more and more severe at higher 
Pt- In fact at large pt the time available for the pair to organize itself into a state with 
the right quantum numbers becomes shorter, and we pretend that it holds together, with 
nothing happening to it which could change its state, until the exchange of Coulomb 
gluons takes over and binds it. This phenomenon manifests itself with a strong form- 
factor-like suppression of the production at large pr, typical of such exclusive processes. 
Since the probability that the quark pair can be found in the right state depends directly 
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on the details of the hard process which produced the pair, it should come as no surprise 
that no universal factorization applies in this case, and that the pt slopes of J/ip and x 
differ. 



It is possible [jrjj to incorporate within the CSM the effect of longer time scales by 
considering higher-order contributions in perturbation theory (PT). For example, one 
could consider production of the heavy quark pair from gluon splitting at large px (as in 
the case of the CEM), and describe the perturbative evolution of this pair into a colour 
singlet state with the right quantum numbers via gluon emission. The gluon virtuality 
before its splitting is of the order of the quarkonium mass, and the time available for 
the quark pair to evolve into the right state is larger than for the LO process. Since we 
allow for the emission of gluons after the creation of the pair, the process is no longer 
exclusive, but rather an inclusive fragmentation process, and the form factor suppression 
is avoided. In the case of Jfip production, diagrams which contribute to the production 
via fragmentation first appear at 0{a 5 s ) in PT (see fig. |TJc). They are therefore suppressed 
by a factor a 2 relative to the LO contributions. However, their p? behaviour is governed 
by the exchange of the gluon in the t-channel, and is therefore given by 1/py. The ratio of 
fragmentation over LO cross sections is then of order (cKsPtA^qq) 2 - This becomes larger 
than 1 as soon as px is larger than few times the quarkonium mass, namely in the region 
where the Tevatron data come from. 

The lesson to be learned is that in the case of quarkonium production, naive a s 
power counting does not establish by itself the correct perturbative expansion, in spite 
of the smallness of a s . In fact, it turns out that the leading-order terms in a s represent 
contributions that should be considered as higher-twist corrections. 

The inclusion of fragmentation contributions bridges the gap between the two extreme 
philosophies of the CEM and of the CSM: in the fragmentation approach, factorization is 
achieved via the separation between the production of a hard (but almost on-shell, relative 
to the global Q 2 of the hard process) gluon and the evolution of this gluon into a given 
quarkonium state. This evolution is described by universal, although state-dependent, 
fragmentation functions, which replace the simple-minded overall constant Ah introduced 
in the CEM. 

Fragmentation functions for all states of interest have been calculated over the past 
couple of years |TJ], [TJ| , and have been used for phenomenological studies . In the case 



of x production, the theoretical calculations agree with the available data, as shown in 
fig. |2| [[17]]. In the case of the 3 Si states, however, the discrepancy in overall normalization 
is striking, although the px distribution fits the data well (see dotted and dashed curves 
in fig. |H). The conclusion is that while fragmentation contributions are fundamental to 
produce the right pr dependence, and are sufficient to correctly predict the absolute 
normalization of the x cross section, there must be additional contributions to explain 
the abundance of ip' produced. 
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Figure 2: Inclusive pt distribution of ?/>'s from Xc production and decay. Upper 
curves and data points correspond to the prompt component. Lower ones corre- 
spond to the b decay contribution. CDF data versus theory. The b contributions 
were evaluated using the NLO matrix elements inclusive b production [18] and 
the CLEO measurements of BR(b — > y) pi 



2.3 The Colour Octet Mechanism 

In order to understand the possible origin of these residual discrepancies, one has to look 
more closely into the structure of the fragmentation process within the CSM. In the case 
of fragmentation into a 3 Si state, the transition of a gluon into the J PC = l + ~ state and 
only one gluon is forbidden. Therefore emission of at least two gluons is required (fig. HJc). 
It turns out from the explicit calculation of the fragmentation probability that only hard 
gluons contribute to this process: soft gluons occupy a small volume of phase space, 
and there is no dynamical enhancement in their emission. Therefore the fragmentation 
probability is proportional to cvf (m). On the contrary, the fragmentation of a gluon into a 
X state requires emission of just one gluon, as allowed by the different quantum numbers 
of the 3 Pj states. Furthermore, emission of a soft gluon is enhanced by the presence of 
a well known logarithmic infrared singularity [pO], which can be regulated by noticing 



that the heavy quarks inside the bound state are slightly off-shell, therefore cutting off 
the singularity at an energy of the order of few hundred MeV. So the probability for a 
gluon to evolve into a x state is of the order of a s (m), as the large logarithm compensates 
the additional power of a s . As a net result, production of x's is significantly enhanced 
relative to that of J/ip's. 
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Figure 3: Inclusive prompt tp' pt distribution. CDF data versus theory. We 
show the contribution of the different sources. Dotted lines: LO production in 
the CSM; dashed lines: gluon and charm fragmentation in the CSM; solid line: 
gluon fragmentation in the colour octet mechanism. 



Again one can interpret this phenomenon using time-scale arguments. Since the evo- 
lution of a gluon into a J/ip in the CSM only involves the emission of hard gluons, the 
time scale for the transition is of the order of their energy in the virtual gluon rest frame. 
In the case of the Xi instead, the enhancement of soft gluon emission indicates that the 
lifetime of the heavy quark pair resulting from the gluon splitting, before it settles into the 
X state, can be very long. Bodwin, Braaten and Lepage developed recently a framework 



1 21], based on non-relativistic QCD, in which such a long-lived state has a non-zero overlap 
with the X- 111 their formulation, a quarkonium wave function is the sum of contributions 
coming from states in the Fock space in which the heavy quark pair is accompanied by 
long-lived gluons. In the specific case of a Xi one nas: 

\ X j) = 0(1)\QQ[ 3 pP}) + 0(v)\QQ[ 3 s[ 8) }g) + ... (11) 

where the upper indices (1) or (8) refer to the colour state of the pair, and v is the 
velocity of the heavy quark in the bound state. The first term in the expansion corre- 
sponds to the standard non-relativistic limit, in which the quarkonium is made just of 
the quark- ant iquark pair. The second term corresponds to a state in which the pair, in 
a colour octet configuration, is accompanied by a gluon (the angular momentum of the 
quark pair is different in the second state, as the gluon itself carries spin). It is precisely 
this second component of the x state that has non-zero overlap with the long-lived pair 
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coming from gluon splitting. Its presence allows the infrared divergence alluded to above 
to be absorbed via a wave function renormalization, in a way which can be rigorously 
defined and extended to higher orders in PT, and which does not have to rely on an 
arbitrary IR cutoff. The final picture we obtain is therefore as follows: the fragmentation 
function for production of a \ state from gluon evolution consists of two pieces. One is 
of order a 2 , and corresponds to the creation of a heavy quark pair with the emission of 
a perturbative gluon, the quark pair being projected on a short time scale on the 3 .Pj^ 
state. The other is of order a s v 2 , and corresponds to the creation of a quark pair in the 
3 S[ 8 * > state. This state will evolve on a long time scale into the required 3 Pj^ state via 
emission of a non-perturbative gluon. The separation of short time scales from long ones 
is arbitrary, but the result is independent of it, as its redefinition only leads to a change 
in the relative importance of the two contributions. The production of x through this 
second component has been named "colour octet mechanism" (COM). Since it turns out 
that, for charmonium, a s is numerically of the same order as v 2 , the two channels are 
competitive. 



Braaten and Fleming |22| suggested that a similar phenomenon might play a key role 
in J/ip production as well. Colour octet states do in fact appear in the expansion of the 
3 S 1 state 0: 

|V>> = 0(1)\QQ[ 3 S[ 1} )} + 0(v)\QQM%) 

+ 0{v 2 ) ( \QQ[ 3 s[ l %g) + \QQ[ l S^ 8) }g) + \QQ[ 3 D { }%g) )+... (12) 



There is a state in this decomposition which can be accessed by a gluon already at 
order a s , namely \QQ[ 3 S[ 8 ^]gg). Creation of this state via gluon fragmentation will be 
suppressed by a factor v 4 , but this can be compensated by the absence of the two extra 
powers of a s that are required for production of the colour singlet state. So once again 
the two processes could be competitive. When one performs the complete calculation, 
numerical factors appear which leave the COM as by far the dominant one. If we take 



for example the integral of the fragmentation functions [[LJ, p2| , we obtain the following 
ratio of probabilities: 

1(8) 



[Prc.bG? - 3 Si)Uet „ OK rr 2 (O^ 

(i) \ 

3 Si/ 



[Prob(<? - 3 S x )) singlet ~ a 2 (m^) { o^y 1 ) 



where the two objects in the last ratio are defined in [21J and can be shown to be pro- 



portional to the square of the wave functions at the origin for the colour octet and colour 
singlet components of the J/ip state, respectively. The factor n 2 arises as a simple conse- 
quence of the phase space for the two additional gluons present in the final state of the 
colour singlet fragmentation. In order for the octet contribution to be larger than the 
colour singlet one by a factor of 50 (the CDF ip' anomaly), it is therefore sufficient that 

(0l%) ~ J/<<>, (14) 

where we used the fact that numerically a s ~ v 2 . This is exactly of the right order of 
magnitude implied by eq. (O). 
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Figure 4: Inclusive i/i pt distribution. Upper curves and data points correspond 
to prompt ip's, after subtraction of the \c contribution. Lower ones correspond 
to the b decay contribution. CDF data versus theory. 



3 Comparison with the Tevatron and SppS data 

The introduction of the COM provides a potentially important new production channel 
that could explain the CDF ip' anomaly. Unfortunately, aside from the generic feature 
that (O^) should be of order v 4 relative to (O^), we have no precise estimate of its 
real value, although it is expected [^] that lattice calculations could provide it one day. 
What can be done, therefore, is to extract (Ofy directly from the data, fitting the CDF 
measurement to the theoretical pr distribution of J ftp and tjj' |23], [HJ . 

The results of the fits are shown in fig. |3| and |j. As had already been shown by Braaten 
and Fleming, the predicted shapes agree very well in the case of the tp'. Now that data 
are also available |J for the production of prompt J fa's (i.e. J/t/i's not coming from either 
b or x decays), we can reach the same conclusion also for the lowest- lying 3 Si state. Had 
we only included the prediction of the CSM, the disagreement with the J ft/) data would 
have been again of the order of 50. The extracted values for the two new parameters (OB,) 
and {0%t) are given in the figures, and it can be verified that they agree with the crude 
estimate given in eq. ([14]). 

Having determined the values of these unknown parameters, one can use them to pre- 
dict rates in different experimental set-ups. In fig. [|, for example, we show the prediction 
for J ftp production in pp collisions at 630 GeV, compared to the UAl data @. In order 
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Figure 5: Inclusive pt distribution of ^'s at 630 GeV. All sources of ip production 
are here included. UA1 data versus theory. The parameters of the theoretical 
calculation take the values fitted on the Tevatron data. 



to match the experimental analysis, which did not separate between different sources of 
ip production, we included the contribution of b and x decays as well. The theory is 
higher by a factor of 2. It is possible that small-x effects ||24j| , which are expected to be 
responsible for the increase of the b production rate at the Tevatron w.r.t. to the NLO 
calculation by an additional factor of 30% relative to 630 GeV [p5| , play an even more 
important role in the case of charm. This would cause our fit to the Tevatron J/ip and ip' 
data to overestimate the values of the parameters. 

A comparison with fixed-target data, where the pt values accessible are much smaller, 
will require additional work, in order to absorb into the NLO parton densities the collinear 



divergences which appear al low pr [26]. Furthermore, at these low p^s one will need to 



take into account the effects due to the intrinsic Fermi motion of initial-state gluons in 



the proton [12 



Another important test of the COM comes from the study of T production. Here the 
colour octet effects are however expected to be smaller, as v is smaller. Nevertheless the 
complete calculation of T rates within the colour singlet model predicts cross sections 
which are smaller than data || [5| by a factor of 3 in the case of IS and 2S states, and by 
a factor of 9 in the case of the 3S. While the factor of 3 discrepancy between the CSM 
prediction and the data could be explained by the addition of the COM contributions, the 
relative factor of 3 discrepancy between the rates of 3S and IS and 2S states is puzzling. 
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Figure 6: Inclusive px distribution of T(15)'s at the Tevatron. Shown are data 
from CDF and the prediction of the CSM, rescaled by a factor of 3.5. The three 
curves correspond to the inclusion of a kx kick of 0, 2 and 3 GeV. 



It is very tempting to assume that a yet unseen third set of Xb states, x'bi exists below 
threshold, and decays radiatively to the 3S state. Various potential model calculations 
of the bottomonium spectrum support this idea. Simple-minded estimates of the x'l 
contributions show that this process can correct the 3S yield by the required factor of 3. 

The shapes of the px distributions are in principle sensitive to the resummation of 
multiple soft gluon emission from the initial state, as the px values probed by the exper- 
iments are small relative to the mass of the T. Figure |6| shows for example the effect of 
an additional kx kick of the order of 2-3 GeV on the LO distributions, compared to the 
T(1S) CDF data. The comparison with the DO data is similar, both in rate and shape. 

Shortly after this Conference, Cho and Leibovich |23j presented a full calculation of 
the contribution of the COM to the T rates, extrapolating to the bb system the values 
of the non-perturbative parameters obtained from the fits to the charmonium data. To 
properly describe the region px < Mr, where most of the data come from and where the 
use of fragmentation functions is not appropriate, these authors performed the complete 
calculation of the LO Feynman diagrams producing the 3 S[ 8 ^ state. Their results indicate 
good agreement with the data, for the IS and 2S states. In the case of the 3S state, a 
residual factor of 3 discrepancy remains if one only includes the decays of the known Xb 
and x'b states. Their detailed estimate of the effect of production and decay of x'l states 
confirms nevertheless the view that their existence would solve even this residual problem. 
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Although they are hard to detect directly, it is not unlikely that future larger statistics 
accumulated at the Tevatron will allow their unambiguous discovery. 



4 Conclusions 

The field of heavy quarkonium production has enormously benefited from the recent 
measurements performed at the Tevatron Collider. Very rarely in the recent past have 
experimental data been so important in guiding the development of a theory. The initial 
discrepancies by almost two orders of magnitude found between earlier models and the 
data have driven theorists to deepen their understanding of the underlying dynamics, and 
have eventually led to a solid framework within which to operate. The inclusion of the 
colour octet mechanism is now viewed as a necessity for the consistency of the theory, 
rather than as an ad hoc theoretical concoction. While the field is still in its infancy, 
and additional progress must be made before a complete theory is formulated and firmer 
predictions can be made, I believe one can conclude that we are on the right track. 
Several calculations are in progress, which will eventually allow us to test the theory more 
thoroughly. 
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